SAGE 测 光 巡天 数据 处 理 方法 研究 
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摘要 SAGE 是 自主 设计 的 能 够 准确 计算 恒星 大 气 参数 以 及 消光 的 新 测 光 系统 。 对 北 天 除 银 盘 外 共计 约 
12,000 平方 度 的 天 区 开展 了 SAGE 系统 测 光 巡天 ， 计 划 获 取 约 5 亿 颗 恒星 的 高 精度 测 光 数 据 。 单 次 曝光 条 
件 下 1000 完备 星 等 usc ~ 17.3，vscg~16.8 (AB 星 等 ) 。 这 为 研究 银河 系 提供 宝贵 的 测 光 资料 。 介 绍 了 巡 
天 专用 的 数据 处 理 程序 的 研究 和 开发 ， 主 要 研究 了 针对 单 幅 图 像 的 快速 自动 化 处 理 过 程 ， 重 点 介绍 数据 改 
正 、 天 体 测量 校正 、 测 光 和 流量 定 标 过 程 ， 以 及 数据 结果 和 数据 质量 检测 等 。 


关键 字 测 光 ; WK; 测 光 系统 ; 数据 处 理 程序 
中 图 分 类 号 : P111 文献 标识 码 : A 文章 编号 : 1672-7673 (2018) 


E 为 了 能 够 更 高 效 、 更 精确 地 获得 恒星 的 大 气 参数 ， 结 合 多 种 测 光 系统 的 优势 ， 并 加 入 全 新 设计 的 滤 光 
e 片 ， 组 成 了 新 测 光 系统 : SAGE 系统 。 该 系统 由 8 个 滤 光 片 组 成 : Strémgren-Crawford u, SAGE v, SDSS 
O sor b DDO51, Ha wider Ha warow〔 依 次 简称 为 usc, vsaor, g, r, i, DDOSI, Haw, Han) 。8 片 滤 光 片 透 过 
率 曲 线 见 图 1， 中 心 波长 和 带宽 等 见 表 1。SAGE 系统 的 详细 介绍 见 文 [1]， 文 [1] 还 详细 介绍 了 SAGE 巡天 
的 方案 和 科学 目标 等 。SAGE 测 光 系统 可 以 提供 关于 恒星 元 素 丰 度 和 银河 系 演 化 (Stellar Abundance and 
= Galactic Evolution, SAGE) 的 大 量 信息 。 自 2013 年 起 计划 进行 Strómgren-Crawford WK”, Z Je UH SE 
SS) 加 优秀 的 SAGE 测 光 系统 进行 原 巡 天 计划 。 

对 北 天 的 SAGE 系统 测 光 巡 天 计划 覆盖 大 约 12,000 平方 度 的 天 区 ， 高 精度 测 光 (S/N=1000) 完 备 星 等 
履 盖 范围 为 usc ~ 11.0 - 17.5 等 ，vsaGg ~ 10.0 - 16.5 等 ，gri ~ 9.0 - 15.5 等 ， 相 当 于 类 太阳 恒星 的 完备 距 
AX ~ 1 kpc。5o 探测 极限 usc ~ 21.5 等 ，vsacE ~ 21 等 ，gri ~ 19.5 等 ， 到 太阳 邻 域 ~ 6.4 kpe 的 距离 。 本 
巡天 于 2015 年 展开 ， 计 划 在 4-5 年 内 完成 测 光 巡天 观测 以 及 流量 定 标 、 天 测定 标 ， 得 到 一 个 深度 均匀 的 
JER SAGE 测 光 星 表 ， 并 利用 该 星 表 开 展 一 系列 以 银河 系 为 主 的 科学 研究 。 


O 


K 1 SAGE 测 光 系统 定义 
Table 1 The SAGE photometric system 


波段 usc VSAGE g r i DDO51 Hay Hay 
中 心 波长 CAD 3425 3950 4770 6231 7625 5130 6563 6563 
带宽 CAD 300 290 1387 1373 1526 154 29 136 
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kia EHI: 2018-03-05; 修订 日 期 : 2018-04-02 

作者 简介 : 郑 捷 ， 男 ， 博 士 研究 生 . 研究 方向 : 天 文 测 光 巡 天 观测 与 数据 处 理 . Email: jiezheng@nao.cas.cn 

通讯 作者 : 赵刚 ， 男 ， 研 究 员 . 研究 方向 : 银河 系 科学 ， 恒 星 物理 .Email: 


201806.00281 v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


Normalized Transmission of SAGE Photometric System 
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图 1 归 一 化 后 的 SAGE 测 光 系 统 滤 光 片 透 过 率 曲 线 


Fig. 1 The normalized transmission curves of the SAGE filters 


1 巡天 观测 


ARRA tic KG Dec > -5° 的 北 天 天 区 ， 并 且 避 开 了 亮 星 密集 的 银 盘 部 分 ( 银 纬 |b|<10°) ,以 避免 了 过 
多 的 亮 星 以 及 图 像 污染 。 此 外 ， 基 特 峰 和 南山 站 秋冬 季 气 象 条 件 较 好 ， 为 集中 观测 秋冬 季 天 区 ， 以 期 在 有 
限 的 时 间 内 获取 较为 优质 的 观测 数据 ， 避 开 了 12 hr < RA. < 18 hr 的 区 域 。 图 2 展示 了 本 巡天 计划 覆盖 的 
约 12,000 平方 度 天 区 。 

为 了 便于 流量 校正 ， 每 个 观测 天 区 和 相 邻 天 区 之 间 都 留 重 夺 区 域 。 在 制定 巡天 规划 时 ， 以 同一 赤 纬 为 
一 个 条 带 ， 划 分 天 区 。 在 同一 赤 纬 ， 相 邻 天 区 之 间 有 约 20% 的 重 钱 ， 而 不 同 赤 纬 之 间 ， 南 北条 带 之 间 也 有 
2 20% He. Al 3 展示 了 Bok 望远镜 6352 号 观测 天 区 与 相 邻 天 区 之 间 的 重 登 情况 ， 南 山 一 米 望 远 镜 观 
测 天 区 重 琶 情况 类 似 ， 但 是 不 同 望远镜 由 于 视 场 大 小 不 同 ， 天 区 划分 并 不 相同 。 此 外 ， 图 中 仅 展 示 了 视 场 
的 外 轮廓 ，CCD 内 部 的 拼接 颖 未 体现 。 

本 巡天 使 用 了 3 台 望 远 镜 共同 完成 观测 工作 : 美国 亚利桑那 大 学 斯 图 尔 德 天 文 台 位 于 基 特 峰 的 博克 
90 英寸 (292.3 m) 望远镜 (The 90-inch Bok Telescope at Steward Observatory of The University of 
Arizona, USA， 简 称 Bok 望远镜 )， 中 国 科 学 院 新 疆 天 文 台南 山 观测 站 的 1 m 大 视 场 望远镜 简称 南山 一 米 
望远镜 ，NOWT) ， 以 及 乌兹别克 斯 坦 科 学 院 元 鲁 伯 天 文 研究 所 迈 丹 内 克 天 文 台 1 m 蔡司 望远镜 (MAO 
Zeiss-1000 Telescope, Maidanak Astronomical Observatory, Ulugh Beg Astronomical Institute, UBAI, fij ff 
MAO 望远镜 ) 。 


一 


1.1 Bok 望远镜 


Bok 望远镜 一 台 
Ney see, ALA 
本 巡天 使 月 


a 


( 赤 经 方向 宽 约 约 166", aA 
慢 速 读 出 模式 下 ， 全 视 场 读 出 需要 约 37 s, ux 
像 元 数 且 3 


AGE Survey Plan 


图 2 SAGE 巡天 计划 覆盖 天 区 
Fig. 2 Overlap between survey fields 
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图 3 
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这 天 天 区 重合 情况 示意 图 


Fig.3 Overlap between survey fields 


未 完全 按照 比例 绘 


ilo 


径 90 英寸 (4 2.3 m) 的 赤道 式 望远镜 ， 隶 属于 美国 亚利桑那 大 学 斯 图 尔 德 天 文 台 
30"57'46'".5， 西 经 111"36'01".6， 海 拔高 度 2017 m， 典 型 视 宁 度 优 于 15". 
H Bok 望远镜 位 于 主 焦点 的 90Prime 大 视 场 相机 ， 终 端 设备 由 
接 而 成 。 视 场 的 大 小 约 为 1. 08”X1.03”， 每 像 元 大 小 (pixel-scale) 约 0.454". CCD 之 间 有 接 缝 
方向 宽 约 约 54") 。 为 提高 读 出 速度 ， 每 个 CCD 由 4 门 并 行 读 出 


噪声 约 6-10 个 电子 。90Prime 的 布局 如 图 4， 图 


4 块 蓝 敏 背 照 式 4K 4K CCD 


。 在 选用 
中 单位 为 
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图 4 90Prime 拼接 式 相机 布局 示意 图 ， 
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位 为 像 元 ， 并 未 完全 按 比 例 绘制 


Fig. 4 Layout of 90Prime, unit is pixel, not to scale 


由 于 Bok 望远镜 口径 相对 另外 两 台 望 远 镜 较 大 ， 集 光 性 能 强 ， 台 址 大 气 透 明度 较 高 ， 而 且 相 机 在 蓝 端 
效率 较 高 ， 所 以 用 于 蓝 端 的 两 个 罕 波 段 usc 和 vsace 的 观测 。 


1.2 南山 1m 望远镜 


中 国 科学 院 新 疆 天 文 台 的 1 m 大 视 场 望远镜 位 于 乌鲁木齐 市 西南 的 南山 观测 站 ， 北 纬 43*16'45".0, 
东经 87°10'38".3， 海 拔高 度 约 2081 m Bl, ^PJ5g T BEZ7g 2.0"。 该 望远镜 为 地 平 式 ， 主 焦点 上 配备 了 大 
视 场 相机 ， 有 效 视 场 约 1.5°X1.5°, AA 4Kx4K 蓝 敏 背 照 式 CCD， 由 4 个 门 进行 并 行 读 出 。 典 型 的 读 出 
时 间 约 40 s， 读 出 噪声 约 8-10 个 电子 。 本 巡天 使 用 了 该 望远镜 配备 的 SDSS g, r, i 滤 光 片 。 望 远 镜 采 用 文 


[种 开发 的 电 控 系统 软件 进行 控 币 


= 


o 


由 于 SDSS 完成 了 北 天 的 大 部 分 天 区 ， 和 本 巡天 计划 天 区 之 间 有 约 9000 平方 度 的 重合 且 深 度 满足 要 
求 。 所 以 只 需要 对 SDSS i tke ZHAI 3000 平方 度 进行 补充 观测 即 可 。 


1.3 MAO 望远镜 


E OIL BRESLER CAA) MAO 望远镜 目前 正在 进行 技术 改造 ， 改 造 之 后 将 具有 37'X 
37' 的 视 场 。 该 望远镜 为 赤道 式 1 m 反射 式 望远镜 ， 位 于 迈 丹 内 克 观 测 站 ， 北 纬 66°53'47", RA 
38°40'22"， 海 拔高 度 约 2593 m[5]。 计 划 在 MAO 望远镜 上 进行 How 和 How 的 观测 。 


1.4 观测 进度 


本 巡天 自 2015 年 秋季 开始 正式 观测 。 截 止 2017 年 底 ， 各 波段 进展 见 表 2. Bok 望远镜 的 usc 和 vsacE 
观测 预计 将 于 2019 年 完成 ， 而 在 南山 1 m 的 g, r, i 波段 在 2018 年 可 以 观测 完成 。MAO 望远镜 目前 正在 


升级 改造 中 ， 预 计 于 2018 年 可 以 开始 进行 观测 。 


表 2 SAGE 巡天 观测 进 


度 (截止 2017 年 底 ) 
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Table 2 The progress of the SAGE survey (by the end of 2017) 


ChinaXiv 合 作 期 刊 


波段 usc VSAGE g r i 
天 区 数 12364 11376 3359 3354 3353 
覆盖 面积 (deg?) 7913 7280 3359 3354 3353 
完成 百分比 ”(%) 69.7 64.2 79.0 78.8 78.8 


2 图 像 预 处 理 


图 像 预 处 理 主要 对 图 像 本 身 进行 改正 ， 消 除 设备 带 来 的 不 一 致 性 ， 以 及 由 设备 引起 的 部 分 图 像 缺 陷 。 


主要 包括 过 扫描 改正 Coverscan) 、 本 底 改 正 (bias). FAE Cflat-field) 、 串 扰 改正 (crosstalk) 等 


步骤 ， 但 是 对 暗流 (dark current) 不 做 改正 。 


2.1 过 扫描 改正 


在 Bok 望远镜 和 南山 Im 望远镜 拍摄 的 图 像 中 ， 都 提供 了 过 扫描 区 ， 过 扫 


x 表达 了 图 像 读 出 时 的 电 


压 水 平 ， 是 针对 本 幅 数 据 的 本 底 。 图 5 展示 了 Bok 望远镜 的 每 个 门 过 扫描 区 和 图 像 区 之 间 的 布局 关系 ， 单 


位 为 像 元 (未 按 比 例 绘制 》。 南 山 1m 望远镜 的 过 扫描 区 分 布 与 此 类 似 ， 但 是 宽度 为 32 像 元 。 
改正 时 ， 程 序 对 每 一 门 图 像 的 过 扫描 区 逐 行 取 中 值 ， 然 后 对 图 像 区 的 每 一 行 数据 ， 减 去 该 行 对 应 的 过 
扫描 区 中 值 ， 完 成 过 扫描 改正 。 由 于 过 扫描 区 数据 起 伏 较 大 ， 因 此 对 求 出 的 中 值 进行 了 高 斯 平滑 。 平 滑 前 


后 的 过 扫描 区 片段 如 图 6。 


对 于 所 有 的 单 次 曝光 图 像 ， 包 括 本 底 、 平 场 、 科 学 图 像 等 ， 都 是 先进 行 过 扫 


处 理 。 


c 
而 
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图 5 Bok 望远镜 过 扫描 区 分 布 示意 图 ， 未 按 比 例 绘制 〈 像 元 ) 


Fig.5 Layout of image and overscan of Bok, not to scale (pixel) 
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图 6 overscan 片段 (红色) 以 及 平滑 后 的 曲线 〈 蓝 色 ) 
Fig. 6 ADU of overscan region (red) and a smoothed curve (blue) 


2.2 AEE CIE 


A ERAT 
Olradgis 


观测 人 员 在 每 个 晚上 的 观测 前 后 各 拍摄 一 组 本 底 ， 通 常 为 各 10 幅 。 对 每 一 幅 本 底 图 像 进行 过 扫描 改 


正 后 ， 对 每 个 像 元 的 多 次 观测 取 中 值得 到 本 底 的 结构 。 在 后 续 的 平 场 图 像 和 科学 图 像 数据 处 理 时 


了 自身 的 过 扫描 区 改正 之 后 ， 都 会 利用 合并 之 后 的 本 底 结 构 进 行 改正 。 


2.3 暗流 


| ， 在 进行 


Bok 望远镜 和 南山 Im 望远镜 都 采用 液 氮 制冷 的 方式 ，CCD 在 冷却 之 后 ， 暗 流 很 小 。 经 测定 ，Bok d 
远 镜 的 CCD 上 暗流 约 为 7.2e/pixel/hr， 而 南山 1m 望远镜 则 不 超过 2e-/pixel/hr。 由 于 本 巡天 采用 的 曝光 时 


间 都 小 于 60 s， 暗 流 带 来 的 影响 相 比 读 出 噪声 等 可 以 忽略 不 计 ， 因 此 未 进行 暗流 改正 。 
2.4 平 场 改正 


2.4.1 Bok 望远镜 平 场 改正 


在 Bok 望远镜 观测 时 ， 在 每 夜 观测 前 后 都 为 每 个 波段 各 拍摄 一 组 (通常 是 10 "85 [ACP (dome 
flat) ， 并 进行 合并 。 圆 项 平 场 的 优点 是 不 受 天 气 条 件 影响 ， 单 幅 ADU 值 高 (通常 控制 在 20,000 至 
30,000) ， 信 品 比 高 ， 能 够 比较 好 地 反映 像 元 之 间 的 差异 性 。 但 是 圆 项 平 场 的 照明 并 不 能 很 好 地 反映 望 远 


镜 上 的 真实 入 射 情 况 。 另 外 一 方面 ， 利 用 每 晚 实际 观 测 得 到 的 科学 图 像 〈 每 晚 大 约 200 到 400 


lii) 进行 


中 值 合并 ， 得 到 当晚 的 超级 平 场 。 超 级 平声 反映 了 真实 的 望远镜 照明 分 布 ， 例 如 比较 均匀 的 入 射 光 场 ， 和 
观测 时 相同 的 光路 等 。 但 是 由 于 单 幅 科学 图 像 的 天 光 背 景 约 在 20-40ADU〈 随 波段 不 同 ) ， 即 使 将 每 个 晚 
上 的 全 部 图 像 合 并 ， 超 级 平 场 的 ADU 值 较 低 ， 信 噪 比 仍 不 高 ， 不 能 反映 像 元 之 间 的 各 像 元 差异 ， 不 适合 


照明 分 布 的 大 斥 度 结构 ， 也 发 挥 了 圆 顶 平 场 的 高 信 噪 比 优 势 ， 改 正 了 CCD 像 元 之 间 的 效率 差异 。 


2.4.2 南山 1m 望远镜 平 场 改 正 


接 用 于 平 场 改 正 。 因 此 ， 参 考 北京 -亚利桑那 巡天 CBeijing-Arizona Sky Skuvey, BASS) 器 的 数据 处 理 方 
法 ， 利 用 二 维 高 斯 平滑 后 的 超级 平 场 对 圆 项 平 场 进行 大 尺度 的 照明 改正 ， 这 样 既 充分 利用 超级 平 场 改正 了 


南山 Im 望远镜 没有 提供 圆 项 平 场 ， 所 以 在 每 晚 的 黄 寻 和 晨光 时 间 ， 对 每 个 波段 各 拍摄 若干 幅 天 光平 
场 。 天 光平 场 兼 顾 了 较 高 的 信 品 比 ， 以 及 真实 的 照明 情况 ， 是 进行 平 场 改正 的 最 佳 选择 。 但 是 天 光平 场 的 
观测 受到 天 气 影响 大 ， 且 观测 时 间 窗 口 短 ， 观 测 次 数 较 少 ， 容 易 出 现 没 有 适合 的 天 光平 场 的 情况 。 同 时 ， 
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由 于 南山 1m 望远镜 视 场 较 大 ， 天 光平 场 也 存在 照明 不 均匀 的 问题 。 因 此 同样 也 需要 利用 当晚 的 超级 平 场 
进行 照明 改正 。 

在 采用 天 光平 场 进行 平 场 改正 时 ， 如 果 遇 到 当天 晚上 未 拍摄 天 光平 场 ， 或 者 拍摄 不 理想 的 情况 下 ， 通 
常 使 用 时 间 最 近 的 天 光平 场 代替 。 


2.5 串扰 及 其 改正 


在 多 门 读 出 的 相机 中 ， 当 图 像 中 出 现 饱 和 星 像 时 ， 会 发 生 串 执 现 象 。 即 当 一 个 门 的 图 像 中 出 现 饱 和 时 ， 
在 其 他 门 读 出 的 图 像 的 相应 位 置 也 会 出 现 镜像 映射 。 通 常 镜像 信号 和 原始 信号 之 间 的 数 之 比 10-4 量 级 。 在 
源 信 号 同一 CCD 的 其 它 门 的 串扰 (Intra-CCD crosstalk) ， 相 关系 数 较 大 且 为 正 相 关 ， 而 其 它 CCD 上 的 串 
扰 (Inter-CCD corsstalk) 则 相关 系数 较 小 ， 且 大 部 分 为 负 相 关 。 

通过 分 析出 现 串 扰 现象 的 图 像 中 串扰 源 信号 和 镜像 信号 之 间 的 相关 系数 ， 然 后 将 镜像 信号 从 原始 图 像 
中 扣 减 。 同 样 使 用 Bok 望远镜 进行 观测 的 BASS 巡天 [中 ， 也 采用 了 类 似 的 方法 。 

(b) 

图 7 展示 了 Bok 望远镜 观测 图 像 上 一 个 典型 的 串扰 现象 ， 以 及 修正 后 的 数据 。 图 7(a) 为 一 亮 星 在 其 他 门 
图 像 相 应 位 置 引 起 的 串扰 ， 图 7(b) 为 改正 后 的 结果 。 其 中 右 下 角 门 为 串扰 原始 信号 ， 上 面 的 两 个 门 为 其 它 
CCD 中 的 门 ， 下 方 3 个 门 则 为 串扰 源 信 号 所 在 CCD 的 其 它 门 。 


33 } w pr S TT 
de turae area: 
uiu 


图 7 改正 前 后 的 串扰 现象 
(a) 原始 数据 ; (b) 改正 后 数据 ;箭头 所 指向 为 串扰 源 以 及 镜像 信号 


Fig. 7 Crosstalk and correction result. 


(a)original; (b)corrected; source and mirrored signal are pointed. 
mil a = + 
3 天 体 测 量 定 标 


天 体 测 量 定 标的 目的 是 找到 图 像 坐 标 系 和 天 球 坐标 系 之 间 的 转换 关系 ， 在 观测 时 指向 的 基础 上 进行 校 
正 ， 最 终 确 定 探测 的 每 个 源 的 精确 天 球 坐 标 。 


3.4 位 置 定 标 和 图 像 扭曲 〈distortion ) 改正 


本 程序 使 用 由 Emmanuel Bertin 开发 的 成 熟 且 被 广泛 使 用 的 天 文 软件 SCAMP7 进 行 天 体 测量 定 标 。 在 
理 中 采用 两 轮 迭 代 ， 从 而 满足 初始 的 较 大 误差 和 最 终 的 较 高 精度 。 第 1 轮 处 理 时 ， 采 用 较 大 的 匹配 误差 
出 ， 这 样 可 以 适应 初始 FITS 文件 头 中 的 信息 不 够 精确 的 情况 ， 第 2 轮 则 采取 较 小 的 匹配 限制 ， 在 第 1 
初步 校正 的 基础 上 ， 能 够 得 到 较为 精确 的 最 终结 果 。 使 用 SCAMP 软件 所 采用 的 重要 参数 见 表 3。 


sr. 


cme 


X 3 SCAMP 重要 参数 设置 


Table 3 The configuration of major parameters for SCAMP 


参数 第 1 轮 取 值 第 2 轮 取 值 说 明 

MATCH Y Y 是 否 进行 模式 匹配 
MATCH NMAX 0 0 匹配 选 星 上 限 ，0 为 自动 选择 
PIXSCALE_MAXERR 2 1.5 像 元 比例 尺 最 大 误差 倍数 
POSANGLE_MAXERR 5.0 2.0 最 大 方位 角 偏 差 ( 度 ) 
POSITION_MAXERR 10.0 1.0 最 大 位 置 偏 差 ( 角 分 ) 
MATCH_RESOL 0 0 匹配 分 辨 率 ，0 为 自动 
MATCH_FLIPPED N N 是 否 允 许 翻转 坐标 轴 匹 配 
CROSSID_RADIUS 25.0 2.0 匹配 半径 ( 角 秒 ) 
SOLVE_ASTROM Y Y 是 否 计算 天 测 结 
PROJECTION_TYPE SAME SAME 投影 类 型 
DISTORT_DEGREES 3 3 图 像 扭 曲 多 项 式 阶 数 


标准 星 方面 ， 采 用 PPMX (Position and Proper Motions eXtended) 星 表 [8 作为 天 测 参 考 星 表 。PPMX 


星 表 收录 了 全 天 约 1800 万 颗 恒 星 的 坐标 ， 空 间 分 布 均匀 。 毕 标 精 度 高 ， 赤 经 和 赤 纬 方向 误差 约 0.02"; AB 
表 中 恒星 亮度 适中 ，85% 以 上 星 等 为 10.0 mag < V « 15.0 mag， 与 本 巡天 的 观测 深度 接近 ， 可 以 很 好 地 与 
观测 图 像 中 的 源 进行 匹配 。 并 且 PPMX 星 表 体积 较 小 ， 便 于 在 观测 现场 进行 快速 数据 处 理 。 后 续 数 据 处 理 
中 计划 改 用 PPMXL[91 等 星 表 进 行 高 精度 数据 处 理 。 


通常 天 体 测量 定 标 的 目的 是 拟 合 图 像 坐 标 (x, ERRERA, 6) 的 转换 公式 ， 然 后 利用 该 公式 计算 探 
测 到 源 的 天 球 坐标 。 

先 将 图 像 坐标 系 转换 为 以 图 像 中 心 为 原点 的 坐标 系 ， 并 用 (u 内 表示 这 个 坐标 系 下 的 坐标 。 在 不 考虑 图 
像 扭曲 的 情况 下 ， 使 用 FITS 文件 头 中 的 CDi j 字段 Gj BUB 1, 2) 所 表达 的 线性 转换 矩阵 CD， 可 以 将 (u,v) 
转换 为 以 角度 为 单位 的 过 渡 平 面 直 角 坐 标 (&, 刀 ， 随 后 根据 球面 投影 计算 对 应 的 天 球 坐 标 (a, 6)， 目 前 使 用 
的 球面 投影 模式 是 正切 投影 (TAN 模式 ) 。 线 性 转换 公式 : 


B =CDx (D. Co es CD22 


由 于 本 巡天 使 用 的 望远镜 相机 视 场 较 大 〈 大 于 1 平方 度 ) ， 并 且 相 机 都 位 于 主 焦 点 ， 焦 面 实 际 上 为 曲 
面 ， 所 以 需要 对 图 像 进行 扭曲 (distortion〉 改正 。 通 常 的 做 法 是 利用 高 阶 多 项 式 进行 拟 合 。 在 多 种 图 像 扭 
曲 表达 方式 中 ， 本 程序 选择 (Simple Imaging Polynomial, SIP]010 作 为 图 像 扭曲 的 表达 形式 。SIP 在 (1) 式 提 
供 的 线性 转换 的 基础 上 ， 利 用 多 项 式 修 正 函数 人 g 对 (uv) 进行 修正 ， 得 到 |: 


b FUP nio (2) 


vtg(uv) 


在 修正 函数 f 和 g 中 ， 用 4 和 Bo 表示 多 项 式 单 元 uPva 的 系数 ， 可 以 写 出 函数 人 g 的 形式 如 (3] 式 。 
HH, A ORDER 和 B ORDER 表示 在 uy v 两 个 方 同 的 阶 数 ， 通 常 二 者 相同 。 经 过 实践 测试 ， 本 这 天 选择 
A ORDER = B ORDER=3。 
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(3) 


天 体 测量 定 标 的 误差 来 自 多 方面 。(1) 星 像 的 中 心 位 置 判断 需要 通过 拟 合 估 计 ， 会 带 来 误差 ，(2) 使 用 


的 参考 星 表 提供 的 坐标 自身 带 有 误差 ， 本 巡天 选用 的 PPMX 星 表 在 杰 经 、 


杰 纬 方向 上 各 有 约 0.02" 的 误差 ; 


(3) 星 表 提供 的 自行 信息 带 有 误差 (4) 在 进行 匹配 时 ， 参 考 星 和 图 像 星 之 间 有 一 定 的 坐标 差 容 忍 度 。 


所 以 ， 最 终 得 到 的 天 体 坐 标 ， 与 其 真实 值 之 间 存 在 误差 。 


3.2.1 天 体 测 量 定 标 外 部 误差 


通过 对 比 图 像 中 探测 到 的 源 与 参考 星 表 PPMX ER) 之 间 的 差异 ， 将 其 作为 天 体 测 量 定 标的 外 部 误 
差 。 图 8 展示 了 随机 选择 的 一 幅 图 像 的 天 体 测 量 定 标 外 部 误差 分 布 图 。 可 以 看 出 天 体 测 量 的 精确 性 比较 高 ， 


标准 差 约 0.1"， 并 且 在 赤 经 、 赤 纬 方向 上 偏差 很 小 。 


N = 780 
ARA: +0.008"+0. 106" 
ADec: —0.029" + Q.142* 
. . 
. 


H 3. 
i eo $e ^ 
ae 
BOO v» 
i 4 
a % 


-0.4 -0.2 0.0 02 04 
A RA (arcsecond) 


图 8 典型 的 天 体 测量 定 标 外 部 误差 


0% 20% 


Percent 


Fig. 8 External error of the astrometric calibration 


3.2.2 天 体 测量 定 标 内 部 误差 


通过 匹配 相 邻 天 区 重 辣 部 分 的 天 体 ， 可 以 得 到 天 体位 置 定 标 的 内 部 误差 。 除 了 相 邻 天 区 的 重 盖 部 分 ， 
还 有 一 些 天 区 进行 了 多 次 观测 ， 也 提供 了 内 部 定 标 误差 的 数据 来 源 。 此 外 ， 同 一 个 天 区 的 不 同 波段 观测 结 
果 ， 也 可 以 用 来 得 到 内 部 误差 。 图 9 展示 了 对 同一 天 区 的 两 次 观测 的 天 体 测 量 定 标 内 部 误差 ， 图 中 可 以 看 
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出 ， 双 向 误差 分 别 为 GRA 
以 忽略 不 计 。 


Percent 


N-1867  "* 
BRA: +0.014¢£0.145* * 
ADec: 30,002" + 0°166+ * 
E * 4 
-- 


A Dec (arcsecond) 
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= 0.014"+0.145" 和 Opec = -0.002"+0.166"， 弥 散 很 小 ， 并 且 在 两 个 方向 的 偏差 可 


-0.4 -0.2 0.0 02 04 0% 20% 
A RA (arcsecond) Percent 
图 9 典型 的 天 体 测量 定 标 内 部 误差 


Fig. 9 Internal errors of the astrometric calibration 


4 流量 定 标 


在 流量 定 标 上 ， 通 常 有 两 种 方式 。 


行 大 气 消光 改正 ， 得 到 视 星 等 ; 


(1) 使 用 测 光 标准 星 ， 在 测 光 夜 进行 多 次 观测 ， 得 到 不 同 大 气质 


量 下 的 消光 情况 ， 并 拟 合 出 大 气 消 光 曲 线 ， 以 此 计算 不 同 大 气质 量 下 的 大 气 消光 值 ， 从 而 对 每 一 幅 图 像 进 
(2) 利用 已 有 的 巡天 星 表 ， 和 待 测 天 区 之 间 重 合 的 部 分 进行 定 标 。 前 者 


比较 复杂 ， 并 且 要 求 在 测 光 夜 进行 ， 但 是 定 标 精 度 高 ， 并 且 独 立 性 好 ， 后 者 操作 上 简单 很 多 ， 但 是 依赖 于 


41 定 标 方法 


己 知 星 表 的 精度 ， 并 且 在 未 直接 重 辣 的 天 区 ， 通 过 重合 区 传递 进行 定 标 会 损失 一 定 精度 。 


对 于 SDSS g. r, i 波段 ， 采 用 APASS(The AAVSO Photometric All-Sky Survey), SDSS, Pan- 
STARRS1 (Panoramic Survey Telescope and Rapid Response System) [2 等 星 表 作为 流量 定 标 参 考 。 
APASS 提供 了 6 和 干 多 万 颗 恒 星 的 测 光 数据 ， 空 间 分 布 上 覆盖 了 全 天 ; 深度 适中 ，10.0 mag < V < 17.0 mag 


的 恒星 超过 90%，70% 以 上 恒星 的 g. r. i 波段 误差 小 于 0.1， 和 我 们 的 观测 深度 相 比 有 重合 ， 


定 标 。 后 期 将 观测 流量 标准 星 ， 提 高 定 标 精度 。 


适合 进行 


而 对 于 usc、vVsacE、DDO51、Haw、Howm 等 5 个 波段 ， 从 经 过 高 精度 流量 定 标 的 恒星 光谱 库 


CALSPEC03] 和 NGSLGI4 中 选择 了 与 本 巡天 天 区 对 


E 合 或 邻近 的 天 


RE) ， 星 等 约 12 至 15 等 ， 且 非 变星 的 恒星 
的 视 星 等 ， 可 以 作为 流量 定 标 标准 星 使 用 。 


区 中 ， 光 谱 型 覆盖 较 全 面 CA 到 K 型 以 及 白 


的 光谱 ， 和 滤 光 片 透 光 率 曲线 进行 卷 积 ， 可 以 得 到 这 些 恒星 
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在 每 个 观测 夜 都 会 对 标准 星 进行 5 至 8 次 观测 ， 然 后 通过 拟 合 得 到 以 〈4) 式 描述 的 大 气 消 光 
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HAG o 


大 气 消光 曲线 是 仪器 星 等 相对 大 气质 量 的 变化 曲线 。 为 方便 分 析 ， 仪 器 星 等 被 统一 归 算 到 1s 曝光 的 仪器 


标 。 


星 等 。 由 于 在 相同 大 气质 量 下 不 同 波长 处 的 消光 略 有 不 同 ， 还 需要 进行 颜色 改正 。 
EAS ke 是 颜色 改正 系数 ， 而 ko 是 仪器 零点 。 然 后 利用 该 消光 曲线 对 当晚 的 科学 观测 数据 进行 流量 


定义 ka 为 大 气质 量 


Ucali 一 Uinst = kua x airmass + kuc x (Ucati E Veali) + kuo 


Veali — Vinst = kya x airmass + kyc x (Heat m Veali) + kyo 


vsace 波段 。 其 中 usc 消光 | 


图 10 和 图 11 是 2017 年 9 H Bok 望远镜 上 观测 数据 得 到 的 消光 曲线 ， 分 别 为 23 日 usc 波段 和 24 H 


(4) 


| 线 为 0.612x airmass + 0.338 mag， 拟 合 残 差 标准 差 为 0.0033 mag, vsace 消光 


曲线 为 0.443 X airmass + 0.241 mag， 拟 合 残 差 标准 差 为 0.0154 mag。 符 合 本 巡天 的 精度 要 求 。 


后 面 将 在 测 光 夜 


Erg 


行 流量 定 标 工 作 ， 每 晚 观测 20-30 次 标准 星 ， 从 而 获得 更 高 的 定 标 精度 ， 预 计 


5-7 个 测 光 夜 完 成 所 有 天 区 的 定 标 。 


4.2 定 标 误差 


如 图 3， 每 个 观测 天 区 和 每 个 相 邻 天 区 之 间 都 有 约 20% 的 


MJD 58019 Band u 


125 | x MJD 58019 Band u (n=6) 
—— 0,612*x 40.338 + 0.0033 


Airmass 
图 10 大 气 消光 曲线 
Fig. 10 The atmospheric distinction curve 
2017-09-23 Bok uSC 


MJD 58020 Band v 


0904 X MD58020Bandv (n=5) 
一 0443'x +0.241 + 0.0154 


Airmass 
图 11 大 气 消光 曲线 
Fig. 11 The atmospheric distinction curve 
2017-09-24 Bok vSAGE 


面积 重 登 ， 利 用 重 登 区 内 匹配 到 的 源 可 以 进 


行 多 次 观测 的 误差 分 析 ， 此 外 还 随机 选择 了 部 分 天 区 进行 多 次 观测 ， 以 分 析 观 测 质量 。 图 12 和 图 13 展示 


了 某 个 天 区 在 两 次 观测 中 定 标的 差异 。 


图 中 红色 曲线 为 每 个 星 等 差异 的 标准 差 。 
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图 12 uSC 波段 流量 定 标 误差 与 星 等 关系 医 
Fig. 12 The error of uSC flux calibration VS mag 
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图 13 vSAGE 波段 流量 定 标 误差 与 星 等 关系 图 
Fig. 13 The error of vSAGE flux calibration VS mag 


5 天 体 测 光 


5.1 图 像 源 探测 


本 巡天 数据 处 理 程 序 调用 成 熟 且 广 泛 应 用 的 测 光 数据 处 理 软件 Source Extractor (以 下 简称 SE)L5] 对 图 
像 进行 测 光 ， 包 括 源 探测 、 源 定位 、 源 累计 探测 器 计数 等 。 测 量 所 用 的 重要 参数 及 取 值 见 表 4。 相 比 传统 
常用 的 图 像 处 理工 具 包 IRAF (the Image Reduction and Analysis Facility) 5651, SE 调用 方便 ， 速 度 快 ， 便 
于 集成 到 数据 处 理 程 序 中 ， 而 且 SE 能 够 一 次 性 完成 源 探 测 以 及 测 光 ， 不 需要 分 多 步 进行 ， 还 可 以 同时 进 
行 多 种 模式 、 多 种 孔径 测 光 ， 效 率 较 高 。 

SE 提供 了 多 种 测 光 模式 ， 例 如 传统 的 圆 孔 径 测 光 ， 根 据 给 定 的 孔径 进行 测 光 ， 得 到 流量 信息 
(FLUX_APER) ， 并 计算 出 对 应 的 仪器 星 等 (MAG_APER) 。 同 时 ， 根 据 天 光 背 景 信息 ， 以 及 流量 的 泊 松 
统计 误差 计算 该 星 的 流量 误差 (FLUXERR_APER) 与 星 等 误差 (MAGERR_APER) 。 如 果 给 定 了 多 个 孔径 ， 
SE 会 依次 进行 不 同 孔径 测 光 ， 并 且 输 出 一 系列 结果 。 此 外 ，SE 还 提供 了 多 种 不 同 模型 对 星 像 进行 测 光 ， 
常用 的 如 自动 孔径 测 光 (MAG_AUTO ) 、 等 亮度 轮廓 测 光 (MAG_ISOCOR) 、Petrosian 模型 测 光 
(MAG PETRO) 等 。 其 中 Petrosian 模型 主要 用 于 河 外 星系 的 测 光 ， 虽 然 河 外 星系 并 非 本 巡天 的 主要 科学 
目标 ， 但 是 对 于 巡天 中 观测 到 的 河 外 星系 进行 尽 可 能 的 记录 和 处 理 ， 也 是 十 分 有 意义 的 工作 。 

除了 测 光 信息 ，SE 也 提供 了 大 量 关 于 星 像 轮廓 的 信息 ， 如 伸 长 率 (ELONGATION) 、 半 高 全 宽 
(FWHM) 、 轮 廓 椭圆 拟 合 的 长 短 轴 、 倾 斜 角 等 ， 还 提供 了 多 种 不 同方 式 得 到 的 星 像 中 心 坐标 。 

参考 SDSS XX. BASS 巡天 CBeijing-Arizona Sky Survey) 等 ， 本 巡天 也 选择 自动 孔径 测 光 模式 作为 
主要 和 输出。 自动 孔径 模式 具有 较 好 的 适应 性 ， 能 够 自动 根据 星 像 的 半 高 全 宽 进 行 选择 孔径 ， 并 且 采 用 李 
孔径 而 非 圆 孔 径 ， 对 恒星 的 流量 能 够 进行 比较 全 面 、 精 确 的 估算 。 与 此 同时 ， 也 提供 多 种 不 同 测 光 模式 
(包括 不 同 孔径 ) 的 测 光 结果 供用 户 选 择 。 此 外 ， 还 对 孔径 测 光 进行 了 生长 曲线 改正 ， 以 期 取得 精度 较 高 
的 测 光 结 
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表 4 Source Extractor 重要 参数 设置 


Table 4 The configuration of major parameters for Source Extractor 
参数 取 值 说 明 
DETECT_MINAREA 5 探测 目标 时 最 小 面积 要 求 
DETECT_THRESH 4.0 探测 目标 时 最 低 探测 下 限 
ANALYSIS THRESH 4.0 分 析 目 标 时 最 低 探 测 下 限 
PHOT_AUTOPARAMS 2.5, 3.5 MAG AUTO 参数 ，Kron 因子 和 最 小 半径 
SATUR_LEVEL 45000.0 像 元 饱和 上 限 1 
MAG_ZEROPOINT 25.0 仪器 星 等 零点 


二 饱和 上 限 为 进行 曝光 时 间 归 一 化 以 及 增益 归 一 化 之 前 的 数值 。 


5.2 孔径 测 光 以 及 孔径 改正 


本 程序 利用 SE 的 APERTURE 测 光 程序 计算 源 的 圆 孔 径 测 光 流量 。 

对 于 孔径 测 光 ， 和 孔径 大 小 的 选择 是 影响 测 光 质量 的 重要 因素 。 孔 径 增 大 时 ， 源 的 流量 能 够 更 多 地 包括 
在 其 中 ， 因 此 对 流量 的 估计 就 更 加 完备 。 但 是 随 着 孔径 的 增 大 ， 天 光 背 景 噪声 也 更 多 地 被 包含 进去 ， 对 暗 
源 的 信 噪 比 影 响 尤为 明显 。 因 此 需要 根据 源 的 亮度 选择 合适 的 孔径 。 另 外 一 方面 ， 为 了 正确 进行 流量 定 标 ， 
同一 图 像 中 的 源 应 以 相同 的 方式 进行 测 光 。 为 了 提高 暗 源 的 信 噪 比 ， 利 用 图 像 中 扳 立 的 、 信 噪 比 高 的 、 非 
饱和 的 亮 星 ， 求 出 仪器 星 等 相对 孔径 的 变化 曲线 ， 即 孔径 生长 曲线 ， 然 后 利用 该 曲线 对 其 他 源 进行 孔径 校 
正 。 通 过 生长 曲线 进行 孔径 改正 测 光 比 单一 孔径 测 光 能 够 较 好 地 提高 处 理 质 量 [07-18]。 

由 于 焦 面 的 不 同 区 域 的 点 扩散 函数 (Point Spread Function, PSF) 不 同 ， 导 致 孔径 生长 曲线 也 不 同 ， 
为 此 本 程序 对 每 个 门 分 别 进行 校正 。 以 Bok 望远镜 为 例 ， 探 测 器 阵列 一 共 由 16 个 门 进行 读 出 ， 图 14 是 一 
组 典型 的 生长 曲线 。 每 个 区 域 表示 对 应 门 ， 绿 色 标 记 表 示 符 合 条 件 的 亮 星 在 不 同 测 光 孔径 的 星 等 ， 红 色 曲 
线 为 采用 的 改正 曲线 。 


V 


AMP #16 || 
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图 14 用 亮 星 得 到 的 每 个 门 的 孔径 生长 曲线 


Fig. 14 The growth curves of brightest stars in each amplifier 
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根据 Bok 望远镜 在 2017 年 9 月 共 15 个 观测 夜 的 观测 数据 ， 对 视 宁 度 、 测 光 零 点 、 天 光 背 景 亮度 等 
SET T Rit AH o 


6.1 ITE (seeing) 


根据 SE 给 出 的 图 像 中 孤立 的 亮 星 的 半 高 全 宽 Cfull width at half maximum, FWHM) 作为 依据 进行 统 
计 分 析 基 特 峰 观测 站 的 视 宁 度 条 件 ， 约 为 1".5+1".0。 考 虑 到 望远镜 和 相机 本 身 对 星 像 的 展 宽 ， 这 个 数值 略 
高 于 站 点 的 真实 视 宁 度 。 图 15 展示 了 上 述 观 测 图 像 的 视 宁 度 分 布 。 


10% 


mmm Seeing: 1.829" x 0.314" 


6% 


Percent 


4% 
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Seeing (arcsecond) 

图 15 2017 fF 9 月 Bok 望远镜 观测 图 像 视 宁 度 分 布 

Fig. 15 The distribution of seeing at Bok telescope in Sep, 2017 


6.2 测 光 零 点 


测 光 改正 零点 是 仪器 星 等 和 视 星 等 之 间 的 改正 常数 ， 反 映 了 观测 时 的 大 气 等 因素 对 流量 的 影响 。 在 进 
行 流量 定 标 时 ， 为 了 避免 不 同 曝光 时 间 带 来 的 影响 ， 所 有 图 像 被 归 一 化 到 一 秒 曝 光 的 流量 。 图 16 是 2017 
年 9 月 观测 数据 的 测 光 零点 的 分 布 ， 其 中 usc 波段 1612 幅 图 像 ，vsacg 波段 1465 幅 。 
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图 16 测 光 零点 分 布 
Fig. 16 The distribution of zeropoint 
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63 天 光 背 景 亮 度 


天 光 背 景 是 影响 观测 质量 ， 尤 其 是 极限 星 等 的 重要 因素 。 图 17 展示 了 usc 和 vsace 波段 的 天 光 背 景 分 
布 。 由 于 对 夜 天 光 的 流量 使 用 恒星 的 测 光 零点 进行 定 标 ， 而 恒星 的 测 光 零点 还 受 观测 时 大 气质 量 的 影响 ， 
所 以 定 标 得 到 的 夜 天 光 流 量 与 实际 略 有 不 同 。 相 比 基 特 峰之 前 的 天 光 背 景 监 测 09]， 亮 度 基本 相似 。 
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图 17 夜 天 光 背 景 亮度 分 布 
Fig. 17 The distribution of night sky brightness 


7 总 结 和 展望 


本 文 介绍 SAGE 测 光 巡天 的 数据 处 理 过 程 和 处 理 程 序 原 理 ， 并 对 数据 处 理 精 度 进 行 了 分 析 。 本 数据 处 
理 程 序 可 以 对 不 同 望 远 镜 取得 的 数据 进行 统一 的 处 理 ， 并 且 得 到 足够 精度 的 数据 。 目 前 本 数据 处 理 程序 主 
要 针对 单 次 曝光 图 像 进行 处 理 ， 多 次 观测 图 像 堆 车 、 相 互 校正 等 等 ， 还 需要 后 续 继续 完善 。 

SAGE 测 光 系统 ， 是 一 个 自主 设计 的 对 恒星 大 气 参 数 非常 敏感 的 系统 。SAGE 测 光 巡天 自 2015 年 开始 ， 
利用 3 台 望 远 镜 协 作 完成 约 5 亿 颗 恒星 的 8 种 颜色 的 观测 。 通 过 测 光 数据 ， 还 将 得 到 高 空间 分 辩 率 的 可 靠 
的 星际 消光 分 布 。 本 巡天 项 目 对 恒星 物理 、 银 河 系 结构 和 演化 等 研究 将 是 非常 宝贵 的 观测 遗产 ， 同 时 也 将 
河 外 天 体 的 研究 提供 大 量 的 观测 数据 。 
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Research on the Data Reduction of the SAGE Photometric Sky Survey 
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Abstract SAGE photometric system is a self-designed system, which can provide accurate stellar atmosphere 
parameters and extinction. The actual 1006 depths of the single-epoch images are usc ~ 17.3 and vsage ~ 16.8 (AB 
magnitude). We are performing a SAGE photometric sky survey of the northern sky, covering about 12,000 
squared degrees and over 500 million stars. This will be valuable data in research of the Milky Way. This paper 
introduces the research and development of the pipeline for the survey, particularly the data reduction process of a 
single exposure. The procedures include bias and flat correction, astrometric calibration, photometry and flux 
calibration. This paper also introduces the main results and data quality produced by the pipeline. 


Keyword Photometry; Survey; Photometric System; Pipeline 
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